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"t^ I Abstract. The observation of cosmic-ray events above the 

/^ . Greisen-Kuzmin-Zatsepin (GZK) cutoff of Eqzk ~ 5 x 10^^ eV 

I challenges orthodox modeling. We discuss a possible solution 

^^ . which uses standard hot Big Bang cosmology and Standard 

I Model (SM) particle physics augmented only by ~ eV neu- 

y-^ • trino masses as suggested by solar, atmospheric, and terres- 

^ I trial neutrino detection. In this scheme, cosmic ray neutri- 

^O . nos from distant, highest-energy sources annihilate resonantly 

on the relic-neutrino background to produce Z-bosons. The 
highly-boosted {jz ^ 10^^) Z's instantly decay, producing "Z- 
bursts" of highly-coUimated jets of hadrons and photons. The 
Q^ i burst content includes, on average, twenty photons and two nu- 

0^ I cleons with super-GZK energies. We show that the probability 

(— I . for each neutrino with energy within a fraction Tz/Mz of the 

pi i' resonant value Eji — 4 (eV/m^) x 10'^^ eV to annihilate within 

JL . the halo of our galactic supercluster is likely within an order of 

Q I magnitude of 1%. Depending on the magnitude of the cosmic 

r*^ . neutrino flux above 10^" eV, this "local" rate for primary pro- 

^ I duction may be high enough to produce the cosmic ray events 

observed above the GZK cutoff. Several tests of this Z-burst 
hypothesis for generating super-GZK events are presented, in- 
^ . eluding (i) a new cutoff energy at En; (ii) a large 7/p ratio for 

" " ' primaries near the upper end of the spectrum; (iii) directional 

pairing of events, and pointing to their cosmic sources; and 
(iv) a neutrino flux above the GZK cutoff energy which is pos- 
sibly measurable directly in proposed ^ 10^^ ton cosmic-ray 
detectors. 
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1. Introduction: the cosmic ray puzzle above 10^" eV 

It has long been anticipated that the highest-energy cosmic-ray pri- 
maries would be protons from outside the galaxy, perhaps produced in 
active galactic nuclei (AGNs). Since the mid-sixties it was also antic- 
ipated that the highest energies for protons arriving at earth would be 
£^GZK ~ 5 X 10^^ eV. The origin of this Greisen-Kuzmin-Zatsepin (GZK) 
cutoff PI is the degradation of the proton energy by the resonant scat- 
tering process p + J2.7K -^ A* -^ N + n when the proton is above the 
resonant threshold for A* production; J2.7K denotes a photon in the 2.7K 
cosmic background radiation. A proton produced at its cosmic source at a 
distance D from earth with an initial energy Ep well above Eqzk will on 
average lose 20% of its energy per interaction length of 6 Mpc to arrive at 
earth with only a fraction ~ {0.8)^^^^^'' of its original energy. Therefore, 
proton energy is not lost significantly only if the highest-energy protons 
come from relatively nearby sources, ~ 50 to 100 Mpc 0. However, no 
AGN sources are known to exist within 100 Mpc of earth.n Any obser- 
vation of air-shower events above 5 x 10^^ eV challenges the theory. To 
date, twenty-four air-shower events with energies near and above 10^" eV 
have been observed by the Volcano Ranch, Haverah Park, Yakutsk, Fly's 
Eye, Akeno/AGASA, and most recently, HiRes collaborations B. It is 
now clear that the highest-energy cosmic-ray spectrum extends beyond 
the GZK cutoff. 

A primary nucleus mitigates the cutoff problem (energy per nucleon 
is reduced by 1/A), but it has additional problems: above ~ 10^^ eV 
nuclei may be photo-dissociated by the 2.7K background, and possibly 
disintegrated by the particle density ambient at the astrophysical source. 
Gamma-rays and neutrinos are other possible primary candidates for the 
highest-energy events. The gamma-ray possibility is not supported by 
the time-development of the Fly's Eye event, but is not ruled out for this 
event. However, the mean free path for a ~ 10^" eV photon to annihilate to 
e+e" on the radio background is believed to be ^ 10 Mpc based on recent 
estimates of the background, and the density profile of the highest-energy 
Yakutsk event showed a large number of muons which may argue against 
gamma-ray initiation. 

1.1. Paired-events and neutrino primaries above 10^° eV 

Turning to the possibility that the primaries may be neutrinos, one en- 
counters an immediate obstacle: the Fly's Eye event occured high in 
the atmosphere, whereas the expected event rate for early development 

^ The suggestion has been made that hot spots of radio galaxies in the supergalactic 
plane (e.g. M87Lat distances of tens of niegaparsecs may be the sources of the super— 
GZK primaries |b|. Statistical support for this hypothesis is weak at present. 
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of a neutrino-induced air-shower is down from that of an electromag- 
netic or hadronic interaction by six orders of magnitude. On the other 
hand, the neutrino hypothesis is supported by evidence [pi that some of 
the highest-energy primaries have common arrival directions, with arrival 
times displaced by only ~ 2 to 3 years. Of the 47 AGASA events above 
4 X 10^^ eV, 9 are contained in three doublets and one triplet with sepa- 
ration angle less than the angular resolution of 2.5°. The chance proba- 
bility of this clustering occuring in an isotropic distribution is less than 
1%. The chance probability for the triplet alone is only 5%. Of the 
seven events above 10^" eV, three are counted among the doublet events. 
Such event-pairing argues for a common source of some duration, emitting 
stable neutral primaries having a small magnetic moment. The stability 
requirement is that ct exceed D ~ Dh, where the Hubble distance is 
Dh = c Hq^ — 1.4 h^^ X 10^* cm, and /igs is the Hubble parameter in 
units of 65 km/s/Mpc. The resulting limit on the rest frame lifetime is 
then trp S> (m/SGZK)(-D_f//c) = 10^^(m/eV)/i^5"^s; for a GeV primary, 
this is about one year, while for an eV neutrino it is as short as tens of 
milliseconds. The requirement that the primary be neutral is necessary to 
allow straight trajectories for the paired primaries as they travel through 
the (probably nanogauss or greater) extragalactic magnetic fields. The re- 
quirement that the primary have a small magnetic moment is necessary to 
avoid significant energy losses via magnetic dipole interactions with the am- 
bient 2.7K microwave background. Neutrino primaries satisfy these three 
criteria. 

Let us call the distance over which a stable particle can propagate with- 
out losing more than an order of magnitude of its energy the GZK distance, 
Dgzk- For a photon it is a few Mpc at 10^° eV, rising to O(IOO) Mpc at 
10^^ eV, with the exact numbers depending on the strengths of the dif- 
fuse radio background. For a proton it is -Dqzk ~ 50 to 100 Mpc. Here 
we propose that the primary particles which propagate across cosmic dis- 
tances above the GZK cutoff energy are neutrinos, which then annihilate 
with relic neutrinos within the GZK zone {D < Dqzk) to create a "local" 
flux of nucleons and photons above Eqzk 0, HI ■ The annihilation rate de- 
pends upon the relic neutrino background reliably predicted by Big Bang 
cosmology, neutrino clustering in gravitational potentials, and the Stan- 
dard Model (SM) of particle physics augmented with neutrino mass values 
suggested by recent oscillation data. For a sufficient cosmic neutrino flux, 
the hypothesis successfully explains the observed air-showers above Eqzk ■ 



2. Z bursts from resonant neutrino annihilation 

It was noted many years ago l^] that the mean free path (mfp), 
Aj = [n,y. aann{vj + Vj — > Z)]^^ , for a cosmic ray neutrino to annihilate 
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at the Z resonance on a background of nonrelativistic relic antineutrinos 
(and vice versa) having mass rvij and density tIjj. is only slightly larger 
than the Hubble size of the Universe, Dh = 1.4 h^^ x 10^'* cm. This 
means that the annihilation probability per cosmic distance of travel may 
be significant and that absorption spectroscopy is a possible means, in 
principle, for determining neutrino masses S. The energy of the neutrino 
annihilating at the peak of the Z-pole is 

E^] = Ml/2m,^ = 4 [cY/m,.) x lO^ieV . (1) 

At a given resonant energy E^: , only relic neutrinos with the j*^ mass 
rrij may annihilate. (We will sometimes use En generically for resonant 
energy with the understanding that there are really three different resonant 
energies, one for each neutrino mass.) 

The invariant energy-averaged annihilation cross section for the process 
^3 + i^j ^ Z is given by the integral over the Z pole. In the SM this is 

/ds 
j^iannis) = 4TrGF/V2 = 4.2 X 10-32cni2, (2) 

with s the square of the energy in the center of momentum frame. When 
neutrino mixing-angles are introduced in the augmented SM, this result 
continues to hold for each neutrino type j in the flavor or mass basis, 
since the annihilation mechanism is a neutral current process. The energy- 
averaged annihilation cross section (cTann) is the effective cross section for 
all neutrinos within ^SEj^/Eji = Tz/Mz = 3% of their peak annihilation 
energy. We refer to neutrinos with resonant flavor j and with energy in the 
resonant range 0.97 E}^: to 1.03 E^}. as "resonant neutrinos." 

Each resonant neutrino annihilation produces a Z boson which imme- 
diately decays (its lifetime is 3 x 10~^^ s in its rest frame). 70% of these 
decays are hadronic, consisting of a particle burst known to include on 
average about one baryon-antibaryon pair, seventeen charged pions, and 
ten neutral pions |^. The ten tt^'s decay to produce twenty high-energy 
photons. We refer to the end product of this Z production and hadronic 
decay as a "Z-burst." The nucleons (we now mean "nucleons" to include 
the antinucleons as well) and photons in the .Z-burst are candidates for the 
primary particles with energy above the GZK-cutoff g . The two nucle- 
ons may be protons or neutrons; if the neutron decays, it simply produces 
a proton with nearly the same energy as the parent neutron. 

If the Z-burst points in the direction of earth and occurs within the 
GZK distance, then one or more of the photons and nucleons in the burst 
may initiate a super-GZK air-shower at earth. The mean multiplicity in 
Z decay is about 30 0. This dilutes the energy per hadron somewhat 
compared to -E^, but it also provides a larger flux per burst. Shown in 
Figure || is a schematic of the Z-burst mechanism within the GZK zone. 
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Figure 1. Schematic diagram showing the production of a Z-burst re- 
sulting from the resonant annihilation of a cosmic-ray neutrino on a 
relic (anti)neutrino. If the Z-burst occurs within the GZK zone (~ 
50 to 100 Mpc) and is directed towards the earth, then photons and nu- 
cleons with energy above the GZK cutoff may arrive at earth and initiate 
super-GZK air-showers. 



2.1. Neutrino mass and resonant energy 

Note that two crucial elements are required for this mechanism to produce 
super-GZK air-showers: the existence of a neutrino flux at ~ 10^^ eV, 
and the existence of a neutrino mass in the 0.1 to a few eV range. Con- 
cerning possible neutrino masses, the simplest explanation for the anoma- 
lous atmospheric-neutrino flavor-ratio llQ] and its zenith-angle depen- 
dence pl| is neutrino oscillations driven by a mass-squared difference of 
i5to^j„-j ~ 10"'^ to 10~^ eV^ ||l2|, which implies a neutrino mass of at least 
0.03 to 0.1 eV. Also, the recent LSND measurement appears to indicate a 
mass-squared difference of Smf^g^^) ^ 0.2 to 4 eV^ 111], from which one 
deduces a neutrino mass of at least 0.5 to 2 eV. 

Upper bounds on the neutrino masses are available from tritium decay 
experiments and from large-scale structure formation studies. From the 
end point spectrum in tritium decay, one infers the upper bound mp < 4.4 
eV for the weighted rms mass mp = (^ \Uej\'^rn3)^ . In a three neutrino 
universe, the smallness of the neutrino m ass-spli ttings deduced from at- 
mospheric and solar oscillations {dm < \/o"^?tm < 0-1 ^V) leads to the 
universal bound m^- < 4.4 eV for all three neutrino masses [O. In a four 
neutrino universe adopted to accommodate the LSND oscillation signal, 
the analogous bound on each of the four neutrino masses is 
m,^ < [(4.4)2 + Smls^j,]^ ~ 5.4 eV @. 

The nonobservation of neutrinoless double-beta decay provides the fur- 
ther upper bound | ^ ■ U^^rrij] < 0.2 eV for Majorana neutrino masses |15|] . 
This bound is restrictive in particular models where the complex mixing 
elements U^j are known. In general, however, cancellations may occur in 
the sum, in which case the bound is compatible with large neutrino masses. 

According to Big Bang cosmology, the fraction of closure density pro- 
vided by the masses of nonrelativistic neutrinos is 

fi^ = 0.025 ft.gg'^ ^(TTi^^./eV). One sees that eV neutrino masses are re- 
quired if neutrino hot dark matter is to contribute in any significant way 
to the evolution of large-scale structure |1^ . Studies of structure forma- 
tion suggest that fi^ ^ 0.15 |1^, in which case no single neutrino mass 
may exceed ~ 6 eV. Particular studies of cosmic structure formation give 
upper limits to m^ in the 2 to 6 eV range p8[ . This is consistent with, and 
complementary to, the tritium bound just discussed. 

According to eqn. p), an upper bound on the neutrino mass of 4 eV 
implies that no .Z-burst energy will fall below 10^^ eV. This is a comfortable 
factor of 20 or so above the GZK cutoff. From the lower bounds on neutrino 
mass deduced from the SuperK and LSND experiments, one gets upper 
bounds on the 2'-burst energy of 10^'^ and 10^^ eV, respectively. It is 
uncanny that the allowed range for Z-burst energies of 10^^ to 10^'^ eV, 
deduced from cosmology, tritium decay, and oscillation experiments, is just 



right for extending the air-shower spectrum an order of magnitude or two 
above the GZK cutoff! 

2.2. Relic neutrino number density. 

The annihilation/Z~burst rate depends on the rehc-neutrino number- 
density. Hot Big Bang cosmology predicts the mean neutrino density of the 
universe. The density of neutrinos with mass below the neutrino decoupling 
temperature ~ 1 MeV is given by a relativistic Fermi-Dirac distribution, 
characterized by the single temperature parameter T^, ~ 1.95K. T^ is a 
factor of (4/ll)^/'^ less than that of the photon temperature T^ = 2.73K as 
a result of photon reheating from the era of e+e^ -^ 77 annihilation. As it 
is the momentum (not the energy) which red-shifts with the temperature, 
the mean momentum of the neutrinos today is pjy ~ 3 T^ ~ 6 x 10~^ eV. 
Consequently, relic neutrinos with a mass exceeding this value have an en- 
ergy dominated by their rest mass. The resulting mean neutrino number 
density is 

(n,^. ) = (3 C(3)/4^') T3 = 54 cm-^ (3) 

for each light flavor j, with an equal number density for each antineutrino 
flavor. The value of C(3) is 1.202. . . . We make two remarks here. The 
first is that the universe should possess an abundance of antineutrinos as 
well as neutrinos. The reason is that the neutrinos and antineutrinos fell 
out of thermal equilibrium while still relativistic (hence the name "hot 
dark matter" ) . This is in contrast to the electron, proton and neutron, 
for which the nonrelativistic Boltzmann factor exponentially suppressed 
the antiparticle. The second remark is that the predicted relic-neutrino 
density is normalized to the relic photon density, which is measured, via 
the temperature relation T^ = Ti'^^- Consequently, the predicted mean 
density of (n^ . ) = 54 cm~'^ in the absence of a neutrino chemical potential 
must be considered firm. The possibility of a nonzero chemical potential 
and its implication for the relic density and annihilation rate are discussed 



in section 2.6 



2.S. Probability for resonant annihilation to a Z -burst 

In a universe where the neutrinos are nonrelativistic and uniformly dis- 
tributed, the mean annihilation length for neutrinos at their resonant en- 
ergy would be 

A, = ((fTann) («.,))-' = 4.4 X lO^^cm, (4) 

which is 30 /igs times the Hubble distance S . A cosmic ray neutrino ar- 
riving at earth from a cosmically distant source will have traversed ap- 
proximately a Hubble distance of space, so its annihilation probability on 
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the relic-neutrino sea is roughly Dh/\j = 3.0/ig5^% (neglecting cosmic 
expansion). 

For a more careful derivation, we let F^.{E^,x) denote the flux of 
the j^^ neutrino flavor, as would be measured at a distance x from the 
source, with energy within dE of E^. The units of F^.{E^,x) are neu- 
trinos/energy/area/time/solid angle. This flux may be quasi-isotropic 
("diffuse"), as might arise from a sum over cosmically-distant sources 
such as AGNs; or it may be highly directional, perhaps pointing back to 
sources within our supergalactic plane. The production rate of Z's with 
energy within dE of E^, per unit length, area, and solid angle, is there- 
fore dFz{E^,x)/dx — aann{Eu)nuj{x)F^.{Ey,x). Integrating this equa- 
tion over the distance D from the emission site to earth, and integrating 
over neutrino energy then gives the total rate of resonant annihilation, i.e. 
Z-burst production (in the narrow resonance approximation), within the 
distance D: 

f ds 
Fz{D) = SF,^ (D) = ErF,^ {Er, 0) I —[1~ exp{~a,„„(s) 5,(Z?)}] , 

(5) 
where Sj{D) = J„ dx n,^.{x) is the neutrino column density from earth to 
the source at distance D. If <Jarm{s)Sj{D) is small compared to one, then 

Fz (D) « Er F,^ (Er , 0) (aann) S, [D] . (6) 

For S{D) <C l/a'arm(s), the rate for Z~burst generation depends linearly 
on the relic-neutrino column density. Writing S{D) as (riiy ) x _D, we make 
contact with our previous simple estimate, namely Fz{D)/ErF^. {Er, 0) = 
D/\j, with \j given in eqn. (0). Also for S{D) ^ l/a'arm(s), the attenua- 
tion of the neutrino flux over the distance from the source to the GZK zone 
will be small, and we may set D in eqn. (O) equal to the GZK distance to 
get the Z-burst rate within the GZK zone. 

For each 50 Mpc of travel through the mean neutrino density, the prob- 
ability for a neutrino with resonant energy to annihilate to a Z-boson 
is 3.6 X lO^''. Since the branching fraction for a Z to decay to hadrons 
is 70%, one part in 4000 of the resonant neutrino flux will be converted 
into a .^-burst containing ultrahigh-energy photons and nucleons within 
Dgzk ^ 50 Mpc of earth in this homogeneous, zero chemical potential 
neutrino universe. 

However, we live in a matter-rich portion of the universe. The Super- 
galactic Cluster, the Local Group, and the Galactic Halo each provides a 
"local potential well" which is expected to trap neutrinos and enhance the 
local relic-neutrino density compared to the universal average. Moreover, 
any neutrino-antineutrino asymmetry (parameterized by a chemical po- 
tential) increases the sum of the neutrino plus antineutrino relic densities. 
Therefore, the value of 0.025% for the probability of a resonant neutrino 



creating a Z-burst within the GZK zone, obtained without neutrino clus- 
tering or asymmetry, is the absolute minimum annihilation probability in 
a Big Bang universe. With clustering and/or asymmetry, the probability 
will be larger. 

2.4- Rate enhancement from neutrino clustering. 

Gravitational attraction will cause neutrinos to cluster somewhat in the 
potential wells of baryonic matter. Let us define some fiducial values for 
cluster sizes (Ls) and the mean separation distance of the clusters (Lss)- 
We take Ds = 20 Mpc for the diameter of the Virgo Supercluster, the 
only supercluster well within our GZK zone, and 100 Mpc, typical of the 
distance across a cosmic void, for the supercluster mean separation distance 
Dss- We take Dc = 5 Mpc for the typical diameter of galactic clusters 
and Dec = 50 Mpc as the typical distance between neighboring galactic 
clusters. For our Local Group of ~ 20 galaxies, we take Dcp = 2 Mpc, 
and Dop-Gp = 20 Mpc for the mean separation distance between groups. 
Dg — 50 kpc is a typical diameter of galactic halos (including our own) 
and Dgg = 1 Mpc is a typical distance between neighboring galaxies. 
With these fiducial values for sizes and separation distances, the respective 
density enhancements (Lss/Lg)^ are of order 100, 10^, 10^, and 10^, for 
the Supercluster, Cluster, Local Group, and Galactic halo, assuming that 
neutrino clustering is more or less as efficient as baryonic clustering for 
these scales. 

The geometrical scaling law that governs the increase in annihila- 
tion rate when neutrinos are clustered rather than distributed uniformly 
throughout the universe is {Lss/Ls)^{Ls/Dh) = L^g/LlDn- If the neu- 
trino clustering scale Lg is less than the GZK distance of 50 to 100 Mpc, 
then the Z-burst production probability within the GZK zone also scales 
with this factor, yielding {LIJDh Ll){3.0hg^ %). Including the 70% 
hadronic branching ratio of the Z, one then gets the probabilities 1.3%, 
2.4%, 1.0%, and 0.2% (independent of /les), for Z-burst production by a 
neutrino of relevant flavor at resonant energy traversing the Supercluster, 
Cluster, Local Group, and Galactic halo, respectively. The larger clusters 
all give a robust probability, within a factor of two of each other, and of 
order of a per cent. This is our main result [^, which we repeat for em- 
phasis: the probability for cosmic ray neutrinos at their resonant energy to 
annihilate within the ^ 50 Mpc zone of earth is likely within an order of 
magnitude of 1 %, with the exact value depending on unknown aspects of 
neutrino mixing and relic neutrino clustering. 

If the neutrino cluster is not isotropic with respect to our preferred posi- 
tion at earth, the Z-burst rate will not be isotropic either. The anisotropic 
rate is easily accounted for by generalizing the distance D from the Z-burst 



to earth to be a vector D. The column density integral encountered ear- 
lier then becomes an integral of n^.(x) along the vector D. There is weak 
evidence that the super-GZK events may correlate with the Supergalactic 
plane ||] . Such a correlation would arise naturally in the model presented 
here if the SG plane provided either the super-GZK neutrino flux or the 
potential well in which neutrinos are clustered (or both). 

We now must discuss the efficiency of neutrino clustering. One expects 
the relic neutrinos to be less clustered than the baryons, especially on scales 
as small as the Galactic halo, for several reasons. First of all, neutrinos 
do not dissipate energy as easily as electrically charged protons do. Sec- 
ondly, neutrinos have a much larger Jeans ("free-streaming") length than 
do baryons at the crucial time when galaxies start to grow nonlinearly. 
And thirdly, Pauli blocking presents a significant barrier to clustering of 
light-mass fermions. As a crude estimate of Pauli blocking, one may use 
the zero temperature Fermi gas as a model of the gravitationally bound 
halo neutrinos. Requiring that the Fermi momentum of the neutrinos not 
exceed the virial velocity <t ^ y/MG/L within the cluster, one gets the 
phase space constraint [Q 

n^j5Acm-^ - 103(m^^/eV)3(cr/200kms"i)^ (7) 

The same constraint holds for np- . The virial velocity within our Galaxy is 
a couple hundred km/sec, whereas virial velocities for rich galactic clusters 
are a thousand km/s or more. Thus it appears that Pauli blocking allows 
significant clustering on the Galactic scale only if m„ ^ 1 eV, but allows 
clustering on the larger scales for m^ ~ 0.1 eV. The free-streaming argu- 
ment also favors clustering on the larger scales. This is just as well, for we 
have shown that it is the larger scales of clustering that give the 0(1%) 
probability for annihilation to a Z~burst. 

2.5. Absorption versus emission mass spectroscopy for neutrinos. 

The original resonant annihilation proposal S focussed on absorption spec- 
troscopy, i.e. a measurement of dips in the high energy neutrino spectrum. 
Measurement of the absorption energies would determine the neutrino 
masses via nij — M'^{N)/2E, where {N) ~ 30 is the mean multiplicity 
in the Z-hwcst and E is an observed air-shower energy. High statistics are 
required to measure absorption dips, for an experiment must map out a 
range of the spectrum. On the other hand, in the absence of background 
events, high statistics are not required in emission spectroscopy. In emis- 
sion spectroscopy, the emitted decay products of the absorbed neutrino are 
measured directly. The Z-burst hypothesis is that the primaries which 
initiate the showers observed above -Bqzk are the emission products of the 
neutrino annihilation process. Furthermore, without the bursts, the GZK 

10 



cutoff prohibits any background events. Thus, by collecting events above 
the GZK cutoff, emission spectroscopy is realized on an event by event 
basis. 

A short annihilation nifp (Xj) is advantageous for an absorption spec- 
troscopic experiment, in that it enhances the depth of the dips. However, 
for emission spectroscopy, where the detected particles must arise within 
the GZK zone, a mfp of order of Dh is optimal — a longer mfp reduces 
the annihilation rate whereas a shorter mfp reduces exponentially the neu- 
trino flux arriving in the GZK zone. Quantitatively, the probability for 
a Z-burst per resonant neutrino is the product of the probability for a 
neutrino to propagate across a distance ~ Dh , times the probability to an- 
nihilate within the GZK distance, i.e. P(Z~-burst)=(e~^«/^0(-^f2Ji), for 

Dgzk <C Aj. This burst probability may be written {^^^){re^^), where 
r = ■^. This probability is maximized at ^ ^^zk -^j^gj^ j. — 2, i.e. when 
Xj = Dh- This derivation assumes a static universe. The result is altered 
by the expansion of the universe. 

2.6. Rate enhancement from a neutrino-antineutrino asymmetry. 

A further enhancement of the neutrino density arises if the universe pos- 
sesses a net neutrino number. Baryon number and lepton number are not 
conserved quantities in the SM, due to coherent processes at the electroweak 
symmetry-breaking scale. Baryon number nonconservation, along with 
CP-violation and out of thermal equilibrium dynamics, may explain the 
observed baryon asymmetry in the universe, AB = [tib — nH)/n^ ^ 10~^. 
Since the baryons are primarily protons, charge neutrality requires a similar 
asymmetry for the electron-positron system. The neutrino asymmetry is 
unknown, and poorly bounded by observation. In terms of the degeneracy 
parameter S, = m/^^? the neutrino asymmetry is 

Az/ = (n,^ - n^^)ln^ = 0.025(71^^ + ^^). (8) 

With a large degeneracy parameter, the thermal equilibrium distributions 
are such that the density of one species {rii, . or np . ) is suppressed ex- 
ponentially in ^ while the density of the other is enhanced as a power- 
law in ^. Also, with C 7^ 0, the sum n,y. + rtp. is always enhanced rela- 
tive to the symmetric ^ = sum. Studies of galaxy formation and pri- 
mordial nucleosynthesis lead to the bounds pOl —0.06 ~ ^^^ ~ 1.1 and 
1^1/ ^ I ~ 6.9. These bounds on ^ translate into Av^ ~ 0.3, which is modest, 
and Ai/^ T ~ 10, which is large. The latter bound allows an enhancement 
in ri^^^ or »T.p^_^ by a factor up to 80. Of course, ^ cannot be so large 
that fiy exceeds the bound deduced from large-scale structure formation. 
Since £,^ {m^/eW) — 430 Q.1, h^^ is valid for m^/Ti, ^ ^ ^ 1, we have, from 
n^ his ^ 0.15, that ^3 (m^/eV) - 65. 
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It has been noted in ||2l| that a degeneracy parameter of ^ ^ 5 (giving a 
neutrino asymmetry factor of '^ 4) effects an enhancement in the neutrino 
or antineutrino number density of ~ 30, thereby decreasing the resonant 
annihilation mfp for antineutrinos or neutrinos to ^ Dh and maximizing 
the Z-burst probabihty within the GZK zone to i x -j-^^ -- 0.2/i65%. 
With neutrino clustering in addition to the asymmetry discussed here, this 
probability is larger still. The factor of ^ in the probability reflects the 
fact that with a large asymmetry, one of n,y or rip is enhanced but the 
other is driven near zero. Consequently, the annihilation probability for 
the cosmic-ray D or v, respectively, is negligible. 

If there is a net neutrino number enhancement of ~ 30, then the 
cosmological bound fli, ^^ 0.15 requires that m^ ~ 0.4 eV for the en- 
hanced species. This in turn puts a lower bound on the resonant energy of 
Er ^ 1022 eV. 

We remark that although the two spin states of Majorana neutrinos are 
not particle and antiparticle, the same thermodynamic counting that ap- 
plies to Dirac neutrinos applies also to Majorana neutrinos. Consequently, 
the asymmetry and density enhancement that result from the calculation 
apply to both kinds of neutrino. The calculation is common because the 
two Majorana spin states can annihilate only with each other when rel- 
ativistic, like particle and antiparticle, owing to the helicity conserving 
nature of the vector and axial vector weak interaction in the relativistic 
limit. Once a density asymmetry or enhancement is fixed at the time of 
decoupling, when the neutrinos are relativistic, it cannot be changed (in 
comoving coordinates) during free expansion. 

2. 7. Two subtleties 

We now discuss two subtleties that will affect the annihilation rate. The 
first is the depolarization of the relic neutrinos as they evolve from their rel- 
ativistic state at decoupling to the nonrelativistic state which they occupy 
today. The second is flavor-mixing among massive neutrinos. 

2.7.1. Neutrino depolarization. The red-shifting of the neutrino momenta 
due to the expanding universe, from relativistic at decoupling to nonrela- 
tivistic today, renders the neutrinos unpolarized. As a result, if the neutrino 
is a Dirac particle, then the densities of the sterile right-handed neutrino 
and the sterile left-handed antineutrino fields are equal to the densities of 
the two active fields. Therefore, for Dirac neutrinos the active densities 
available for annihilation with the incident cosmic-ray neutrino are half of 
the total densities, and the Z-burst production probability is half of what 
we quote in this article. For Majorana neutrinos, there are no sterile fields 
and the total densities are active in annihilation. 
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Depolarization of Majorana neutrinos also means that rt^^ — 7i,y^ today. 
Although the concept of "neutrino asymmetry" or "net lepton number" 
loses its meaning for Majorana neutrinos, any density enhancement due to 
an asymmetry at the time of decoupling remains. Majorana neutrinos are 
favored over Dirac neutrinos in currently popular theoretical models with 
nonzero neutrino mass 123]. 

2.7.2. Flavor mixing in the mass basis. Oscillation experiments strongly 
suggest that the mass and interaction ( "flavor" ) bases of neutrinos do not 
coincide. In general, the flavor states are unitary mixtures of the mass 
states. Letting a = e,/u,T label flavor and j = 1,2,3 label mass, one has 
\i'a >= ^ Uajl^j >, where U is the unitary mixing matrix. Then each 
neutrino flavor at the resonant energy of a given mass state has a nonzero 
probability to annihilate, but with an extra probability factor of |C/ajp. 
For example, the f^'s and I'e's from pion and mu decay will annihilate at 
the resonant energy of m2 with the probability factors | C/^2 1 ^ and | Ue2 P , re- 
spectively, times what we have calculated above. The mixing factors |C/ajP 
can be easily multiplied into our result. Indications from the solar and at- 
mospheric neutrino experiments are that the mixing factors are large; e.g., 
in the presently popular "bi-maximal" model [g3|, all the |C/aj"P factors, 
except for |C/e3p) are in the range 0.25 to 0.5. 

Flavor mixing also affects the present density of relic mass eigenstates, 
giving n^^ — J2a Wajl'^'n^^, where J2a ^^ ^ ^^^ '^^^^ ^^^ three neutrino 
flavors and n^^ is the Va relic density expected in the absence of mixing. 
This formula results from the separation ( "decoherence" ) of the neutrino 
wave function into distinct wave packets for each mass eigenstate after a 
Hubble-time of travel. Note that if the three Va are the same, = n,y, then 
unitarity gives n^. — n^, for all three mass states j — 1,2,3. 

Finally, we mention that the phenomenon of neutrino oscillations oc- 
curs when neutrinos are mixed. Neutrino oscillations will affect the flavor 
populations of the cosmic-ray neutrinos. However, they will not affect our 
calculation of annihilation. 



3. Cosmic neutrino flux above 10^" eV 

A considerable cosmic neutrino flux above Eqzk is required for the Z- 
burst hypothesis to successfully explain the super-GZK events. According 
to eqn. (^), the requirement on the neutrino flux at the resonant energy is 
that the product of this flux per flavor times the resonant energy, times the 
annihilation probability within the GZK zone (which we have estimated to 



be 10 ^^^ in section 2.4), times the photon and nucleon multiplicity per 



burst (^ 30), is comparable to the flux of events above Eqzk', this is, 

ErF,^{Er) ^ 10°-^±iFp/^(> ^gzk). (9) 
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The AGASA collaboration |2J shows an integrated flux above 5 x 10^^ eV 
of Fp/^(> 5 X IQi^eV) - 1.0 x IQ-i^/cniVs/sr. Based on the twenty 
events above 10^° eV with the best energy calibration,^ A. Watson p5[ 
calculates an integrated flux above 10^" eV of Fp/^{> 10^° eV) = 2.0 ± 
0.5 X 10~^°/cm^/s/sr. These two fluxes are quite consistent with the spec- 
tral index of 2.8 ± 0.3 found for energies above 10^^ eV by the AGASA 
collaboration [e3 . It seems safe to assume that the flux Fp/^ which is gen- 
erated by the Z-bursts must be bracketed by these two flux values, i.e., 

by 

Fp/,{> Egzk) = 10-i9-5±°VcmVs/sr. (10) 

This allows for some of the events just above Eqzk to arise from relatively 
nearby sources via a conventional mechanism, or from statistical fluctu- 
ations or systematic errors. From eqns. (H) and ([lO|), we arrive at the 
following estimate for the resonant neutrino flux: 

EB.F,,{En) ^ 10-i9±iVcmVs/sr. (11) 

It is not unlikely that whatever mechanism produces the most energetic 
hadrons also produces charged pions of comparable energy. Thus, one may 
expect neutrino production at ultrahigh energy, coming from pion decay 
and subsequent muon decay. Photopion production in a low-density source 
with a straightforward power-law extrapolation of the proton spectrum is 
not promising p6[ . However, a dense source such as an AGN may effec- 
tively trap hadrons so that the confined proton flux greatly exceeds a power 
law extrapolated from the observed cosmic ray flux p^ . The possibility of 
a large high-energy neutrino flux does not seem unnatural. In fact, it has 
been claimed that the MHD "snowplow" effect in a dense plasma is capable 
of producing an enhanced neutrino flux all the way up to 10^"* eV, with- 
out the concomitant proton and photon fluxes violating any observational 
bounds m. 

There is also the possibility that the highest-energy neutrinos originate 
in quark jets produced in the decay of some supermassive relic particles p9| 
or topological defects p^| , in which case the neutrino flux greatly exceeds 
the proton flux. However, this kind of exotic source takes us beyond SM 
physics, thereby mitigating the economy of the Z-burst hypothesis. 

There do exist experimental upper limits on the high energy neutrino 
flux, derived from the nonobservation by the Fly's Eye of air-showers pro- 
duced by penetrating primaries Q. Using the QCD-corrected high-energy 
neutrino cross-section estimated to be ct^at ~ 8 x 10~'^^£^26* '^^ [eH' with 
E20 being the neutrino energy in units of 10^° eV, I get the following for 

^ Watson's flux is based on the twenty events above 10^" eV from Volcano Ranch 
(1), Haverah Park (4), Fly's Eye (1), AGASA (7), and HiRes (7); the Yakutsk events 
are not included here since there is some uncertainty in their energy calibration. 
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the Fly's Eye bounds: 2 x IQ-^^, 3 x 10"^'', and 5 x 10"^^ per cm^-s-sr 
at Ely = 10-^"^, 10"^^, and 10^'^ eV, respectively. Thus, the neutrmo flux at 
10"'^^ eV cannot be much more than six orders of magnitude larger than the 
flux at the resonant energy. At present, this is a constraint on the source, 
not on the Z-burst model. 



4. The p, n, 7 flux above Eqzk from Z bursts 

The decay products of the Z are well-known from the millions of Z's pro- 
duced at LEP and at the SLC 0]. Among the ~ 30 particles in a Z-burst 
are 20 photons from 7r° decays, and 2 nucleons; these are the candidate 
primary particles for inducing super-GZK air-showers in the earth's atmo- 
sphere. 

Without simulation and modeling it is not clear what ratio of photon 
to nucleon initiated air-showers is to be expected above the GZK cutoff. 
The a priori photon-to-nucleon ratio is about 10. However, the hardness 
of the nucleon spectrum compared to the photon spectrum as measured 
in Z-decay mitigates this ratio if a selection is made for the very highest 
energy particles. The mean pion energy from Z decay is about a third of 
the mean nucleon energy M. The mean photon energy from 7r°-decay is 
half again that of the pion. Complicating the predicition for the photon 
to proton ratio for air-shower initiation is the different physics causing the 
attenuation of photons on the radio background, and degradation of the 
proton energy on the 2.7K microwave background. Also, the development 
of photon-initiated air-showers above Eqzk is skewed by the coherent elec- 
tromagnetic "LPM" effect and by high-altitude photon-absorption on the 
earth's magnetic field pS, either of which may affect the shower identifi- 
cation and interpretation. Finally, the average values for the multiplicities 
and energies presented here must be used with some caution, since fluc- 
tuations in multiplicity and particle- types per event, and in energy per 
individual particle, are large. For example, although the mean number of 
baryon-antibaryon pairs per hadronic Z-decay is one, the probabilities for 
0, 1, 2, 3, and 4 pairs are 37%, 37%, 18%, 6%, and 1.5%, respectively, if 
the number of pairs is governed by Poisson statistics. 

What does seem clear is that there should be a dominant fraction of 
photon-initiated air-showers not far below the upper end of the observed 
spectrum. Two recent numerical simulations of the Z-burst hypothesis 
have been performed, including updated quark-to-hadron fragmentation 
functions at the Z-resonance and all known extragalactic propagation ef- 
fects. The conclusion of the first study [Q is that if the high energy 
neutrinos are emitted from a dense source with a spectrum extending to 
the resonant energy, and if relic neutrinos of ^ eV mass are clustered on 
the Supercluster scale by a sensible amount, then the Z-hwcst mechanism 
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can produce the events and rate observed above the GZK cutoff without 
violating lower energy ffux limits. A further simulation result is that the 
fractions of photon and nucleon primaries are similar above 10^*^ eV, with 



photons dominating below 10^° eV. The second study |35 admits the pos- 
sibility of a successful Z-burst explanation of the super-GZK events as 
long as extrapolation of the cosmic neutrino flux down in energy satisfies 
the Fly's Eye bounds discussed in section 0. 



5. Further signatures from Z bursts 

The Z-burst hypothesis is highly predictive and testable. Here we discuss 
some possible signatures: 

(i) As just mentioned in section U, the 7/p ratio for air-shower primaries 
should be large not too far below the end of the spectrum, 
(ii) The energy of the Z-bursts are fixed at 4 x [eW /m^.) x 10^^ eV by the 
neutrino mass(es). The energy of individual particles produced in the burst 
can approach this value but cannot exceed it. This may serve to distin- 
guish the Z-burst hypothesis from some recent speculations for super-GZK 
events based on SUSY or GUT-scale physics, in which cutoff energies are 
expected to be much higher. 

(iii) From the highest super-GZK event energy E"'^^^, one can deduce an 
upper bound on the neutrino mass of 

m^ < Af|/2£:"^''^ = 4 {lO^^eV/E"''''') eV . (12) 

Similarly, from the mean energy (E) of super-GZK events one can estimate 
the mass of the participating neutrino flavor via 

m^ = MI/2Er - Ml/{2{N){E)) - 1.3 x {IQ^° eV / {E)) eV . (13) 

If there is a selection bias toward events at higher energy, then this formula 
gives a lower bound on the neutrino mass. 

(iv) The highest-energy neutrino cosmic-ray flux should point back to its 
sources of origin, and the super-GZK event arrival directions should point 
back to these same sources. The AGASA observation of directional pairing 
(section |l.lD of highest-energy events suggests that this is happening. In- 
deed, a possible correlation between arrival directions of the E ^ 10^° eV 
events and the locations of compact radio quasars has recently been noted 

The gamma primaries from Z-bursts must point straight back to their 
source, whereas the nucleon primaries from Z-bursts will be somewhat 
bent by the magnetic fields in the GZK zone. After a randow walk through 
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a distance D of magnetic domains with coherence length A, the nucleon 
bends through angle 

50 ^ 0.5° X (Dmpc Ampc) ^ B„G E20 . (14) 

Here B^g is the magnetic field in nanogauss. So, 10^" eV nucleons coming 
from 50 Mpc will bend a few degrees if they encounter nanogauss domains 
of Mpc size. Conventional wisdom is that extragalactic fields of strength 
nanogauss or less are to be expected [M. However, there is a recent claim 
that the extragalactic fields may actually be of microgauss strength |3q] . If 
true, then nucleons produced in Z-bursts will bend dramatically, and bear 
no directional correlation whatsoever with the cosmic sources. 

Note that neutrons will bend as much as protons for two reasons. First, 
for each 6 Mpc (mfp for photopion production off the 2.7K background) 
of travel, the proton and neutron have an equal chance to become each 
other via charged pion emission. Secondly, for each neutron decay length 
CTn = E20 Mpc traveled, the neutron has a probability 1 — e~^ = 63% 
to decay into a proton with negligible energy lost to the accompanying 
electron and neutrino. 

(v) There could be a "neutrino pile-up" two decades of energy below En. 
These pile-up neutrinos are the result of the hadronic decay chain which 
includes Z ^ ~ 17 {tt^ — > f;j + /i^ ^ t'/x + f'/j + Ve/i^e + e^); *-e. each of 
the 70% of the resonant neutrino interactions which yield hadrons produces 
about 50 neutrinos with mean energy ^ -E7r/4 ^ -Ei^/120. These neutrinos 
are in addition to the neutrinos piling up from the decay of pions photo- 
produced by any super-GZK nucleons scattering on the 2.7K background. 
(vi) The Lorentz factor of the bursting Z is 

7z = Er/Mz = Mz/2m, = 0.9 {m,/eV)-' x lO". (15) 

The Z-decay products which in the Z rest frame lie within the forward 
hemisphere are boosted into a highly-coUimated lab-frame cone of half- 
angle 1/jz = 2(TO^./eV) X 10""'^^ radians. Z-bursts originating within 
20 {eV/m,y- ) parsecs of earth, if directed toward the earth, arrive with a 
transverse spatial spread of less than one earth diameter. It is therefore 
possible (but unlikely) for the decay products of a single Z-burst to ini- 
tiate multiple air-showers. A large area surface array {e.g. the Auger or 
Telescope Array projects) or an orbiting all-earth observing satellite {e.g. 
the OWL/AirWatch proposal) could search for these coincident showers. 

There are a few more remarks of interest concerning the Z-burst hy- 
pothesis: 

(i) If the Z-burst hypothesis is correct, then the neutrino flux at _E ~ 10^^ 
eV is sufficiently high that a direct measurement of it is possible with a 
teraton (lO^^ ^n) detector. With a^N ~ 8 x lO'^^ E^q^ crn^ [||, and the 
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neutrino flux given in eqn. (pi]), one has (1 ton = 0.6 x 10'^" nucleons): 

Events/year/sr/ton = 10-i3-2±i.4 (^^)"-i J dE,, E'i,'-''^ (16) 

within the energy range given by the integration hmits. Here n is the 
spectral index of the neutrino flux in the vicinity of the resonant energy. 
Invoking as an example a flat spectrum (n — 0) above -Eqzk up to i^maxj 
one gets ~ lO""^^^'* (%f^)(i§S^)°'' events/yr/sr/ton above the GZK 
cutoff. For the possible values i?max ~ 10 Efi ^ 10^'^ eV, the rate from this 
flat spectrum example is ~ ^.O^^"'**^'* events/yr/sr/ton, which is probably 
observable in a teraton detector. Whereas the volumes of operational de- 
tectors are ~ 10® tons for AGASA and 10^" tons for HiRes, the volumes for 
the larger of the proposed detectors are 10^^ tons for Auger and the Tele- 
scope Array, and 10^^""'^^ tons for OWL/AirWatch. Since neutrinos with 
energy above a PeV have a charged-current interaction length in matter 
less than the earth's diameter |39|, the neutrino signature of interest in 
these detectors is a penetrating horizontal shower. 

A novel proposal for measuring the highest-energy neutrino flux mon- 
itors for radio pulses from the limb of the moon, which may be produced 
by penetrating high energy neutrinos |Q . The idea is to use radio to "see" 
the interactions of neutrinos traversing a small column-density of matter. 
For detectors in earth orbit, neutrino interactions in the limb of the earth 
could provide an analogous radio signal. 

(ii) If the highest-energy neutrino flux is nearly isotropic, then the super- 
GZK event directions should correlate with the spatial distribution of the 
relic neutrino density. The solid angles subtended by any nearby halos may 



offer preferred directions for super-GZK events. As discussed is section 2.4, 
the Supergalactic plane may be the most probable cluster domain for neu- 
trinos. Perhaps the angular distribution of super-GZK events can be used 
to perform neutrino-cluster tomography. 

(iii) If the super-GZK events are due to neutrino annihilation on relics 
as hypothesized here, and if the high-energy neutrino flux is eventually 
measured, then an estimate of the relic-neutrino column density out to 
Dqzk ^ 50 to 100 Mpc may be made. Solving eqn. (||) for Sj{D(jzk) and 
replacing Fz{Dqzk) with Fp/^{> Eqzk)/{N), we have 

n fr) \ Fp/i{> -E^GZk) 

Ji(^GZKJ - 



ERF,^{ER){a^^,,){N) 
= O.SxlO-x^kl^^cm-. (17) 

Here {N) is the mean number of nucleons and protons in the Z-burst which 
are above Eqzk- For the numerical example we have taken (N) ~ 30, which 
applies if all the nucleons and photons in the Z-burst are above Eqzk ■ An 
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estimate of Fp/j{> Eqzk) based on recent data was given in cqn. (|1 1|) . 
The feasibility of a future measurement of EjiFu^Ee) was discussed in 
section |[ (If F^j (E) is measured below the resonant energy, an estimate 
of the neutrino column density can still be made by model extrapolation of 
the flux to En.) Thus, an experimental determination of Sj{Dqzk) seems 
possible in the not too distant future. 



6. Summary and prospects 

In summary, if one or more neutrino mass is within the range ~ 0.04 — 
1.0 eV, and if there is a sufficient flux of cosmic ray neutrinos at ^ 10^^ 
eV, then Vci + iridic (or vice versa) ^ Z — > nucleons and photons within 
the GZK volume ~ (50Mpc)'^ of earth may be the origin of air-shower 
events observed above the GZK cutoff. If the hypothesis is correct, then 
air-shower observations may have already shown the existence of the relic- 
neutrino gas liberated from the primordial early-universe plasma when the 
universe was only one second old. 

There are good prospects for more cosmic ray data soon at the highest 
energies. The AGASA and HiRes experiments are active. In the near 
future. Auger and the Telescope Array will become operational and increase 
the sensitivity by a factor of 10 to 100. In the more distant future, the 
proposed OWL/AirWatch satellite experiment, with fluorescence detectors 
looking down at our atmosphere, may become operational and increase the 
sensitivity by yet another factor of 10 to 100. 

Possible signatures to validate or invalidate the Z-burst hypothesis pre- 
sented here will be forthcoming. Signatures are abundant. Several were 
listed in section ||. The most striking of these are: 
(i) a new cutoff energy at Er = 4 (eV/m,y) x 10^^ eV; 
(ii) a large j/p ratio for primaries near the upper end of the observable 
cosmic-ray spectrum; 

(iii) pairing of events in direction, and directional pointing of shower-axes 
to their cosmic sources; 

and (iv) a neutrino flux sufficiently large above the GZK cutoff energy that 
direct measurement is possible with the proposed Auger, Telescope Array, 
and OWL/AirWatch detectors. 
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